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Abstract 

Let T : X — > X be a compact linear (or more generally affine) operator 
from a Banach space into itself. For each x € X, the sequence of iterates 
T n x, n = 0, 1, ■ ■ • and its averages i ]Cfc=o T k ~ 1 x, n = 0, 1, • • • are either 
bounded or approach infinity. 
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Let X be a set and / : X —> X a map from X into X. For an x £ X, 
the sequence of iterations x, f(x), f 2 (x), • ■ • , f k (x), ■ ■ ■ can be considered as a 
trajectory of a dynamical system where time is the discrete nonnegative integers: 
starting with the initial (time t = 0) state x, the state at time t = k is f k (x). 
Suppose now X — C™ and / is the transformation defined by an n x n complex 
matrix A. What can one say about the general behavior of trajectories of A? 
More generally, we consider affine maps on X, i.e. maps of the form Ax + c, 
where A is linear and c is a constant vector, and we also allow X to be infinite 
dimensional. Moreover, we study the behavior of the sequence of averages: 

Ave fc /(x) =Ux + f{x) + ■■■+ f k - l {x)),k = 1, 2, • • • 
k 

Note that the method of averaging was used in [2] in approximating solutions 
of a system of linear equations. In case of linear operators, problems about the 
linear span of the iterates T k x,n = 0,1, •• • can be found in [4]. 
Recall that a square matrix N is called nilpotent if N s = for some nonnegative 
integer s. 

Throughout this paper, for nonnegative integers k,j, k > j, C(k,j) denotes the 
binomial coefficient 

kl 
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By convention, C(k, 0) = 1 for k > and C(k,j) = if k < j. 

If k < j, the sum Y]j—j Ui is considered as an empty sum and its value is 0. 

Theorem 1 Let T : C™ — > C™ 6e on offine mop defined by Tx = Ax + c, where 
A is an nxn complex matrix and c a constant vector in C n . Let || • || be a norm 
on C n . Then for any vector i€C", the sequence 

T k x,k = 0,1,2, ■■■ 

is either bounded or linifc ||T fe a;|| = oo. 

Proof. 

By Jordan canonical decomposition theorem, C™ = V\ © V2 © • • • © V m for some 
subspaccs Vi,i — 1,2, ■ ■ ■ ,m with the following properties: (a) each Vi is an 
invariant subspace of A, i.e. Av € Vi for all v £ Vi, and (b) there exists Aj e C 
and a nilpotent matrix TV^ such that Av = XiV + NiV for all «eKj. 
Let Pj be the algebraic projection of C™ onto Vi associated with the decompo- 
sition C™ = Vi © v 2 © • • • © v m . 

Define a new norm | • | on C™ by 

\ v \ = \\P lV \\ + \\P 2V \\ + ... + \\P mV \\. 

Let a; be a vector in C™. Let 

x k = T k x = A k x + c + Ac-\ h A k ~ 1 c, k = 1, 2, • • • . 

For any vector x, we have the following equalities: 

x = P Y x -\ h P m x 

A k x = A k P lX + ■■■ + A k P m x 
A k x = P x A k x + ■■■ + P m A k x 

The last equality follows from the commutative property of A and Pi , i = 
1, ■ ■ ■ , m since A is invariant in each Vj. 

Fix an i and write v = PiX,d = PiC, X = Xi,V = Vi,N = A/j. Let s be the 
smallest nonnegativc integer such that N s v — and t the smallest nonnegative 
integer such that N l d = 0. Then for k > max{s,t} and s,t > 1 (if s = or 
t = 0, the corresponding sum below is defined as 0), one has 

P iXk = (XL + N) k v + d+(XL + N)d+--- + (XI + N) k - 1 d (1) 

s-l t-1 

= ^C(fc,j)A fe ^'iV^ + ^^,fc)iV^ (2) 
j=o j=o 

where 

S(j, k) = C(j,j) + C(j + l,j)X + ... + C(k- 1, jJA*- 1 ^'. 
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Using the identity C(j, i + 1) + C(j,i) = C(j + l,i + 1), we have for A = 1, 
= C(k,j + 1)) an d f° r 1, fc) is given recursively by 

5(0 ' fc) = T^X 

and (by subtracting XS(j,k) from S(j,k)) 

(1 - A)S(j, fc) = S(j - 1, fc) - X k - j C(k,j),j = 1,2, • • • ,t - 1, ( for t > 2), 
from which we get an alternate formula for fc): 

s( ik )- l - xk y M^r! ... t_i rs) 

6U.«;- (i_ A ),+i 2_, (i-A)J"-< ,J ' ' ( j 

Note that (3) is also valid for j = for ^ s an em Pty sum. We shall show 

that for A 7^ 1, 

IU — 1 

P i x k = B+Y J \ k ~ j C{k,j)A j (4) 
j=o 

where w = max{s, t}, and Aj,B,j = 0, • • • , w — 1 are constant vectors indepen- 
dent of fc; and for A = 1, 

i 

P l x k = v + J2C(k,j)B J (5) 

where I = max{s — l,t}, and Bj,j — 1, ■ ■ ■ , / are constant vectors independent 
of fc. More precisely, for fc > max{s,t}, 

t-i 

^ = - £ +1 JVM,j = 0, • • ■ ,min{M} - 1, 

i=j 

t-1 1 

= e(s - t)N 3 v - e(t - s) ^ — — — - — — J - min{s, t}, ■ ■ ■ , w - 1 



2=J 



and 



t-i 



and 

= N j v + N j - 1 d,j = l,-..,min{s- l,t}, 
Bj = e(s - 1 -t)N j v + e{t- s + l)N j - 1 d,j = min{s - 1, i} + 1, • • • , I, 
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where e(r) = 1 for r > and e(r) = for r < 0. 
Indeed substituting (3) into (2), we get 

i=o j=o \ v ; »=o v ; / 

= E^a-^+e 

j=0 j=o y ' j=0 i=0 v ; 

s 1 t ■ 1 i X , — | / 7 ■ \ \ /c i 

= E i)**- + B - E E a - AV-^i ^ 

j=0 i=0 j=i ^ ' 

= g c(*,i)A*-Wi, +B-EE ^ 

i=o j=o i=j y 1 

s-1 t-1 t-1 

So x 



E C{k,j)X k -iNiv + B-E C(k,j)X k -' E fl _ A y- 3 - + i ^ 



Now (4) follows by considering cases s > t, s = t or s < t. 

Next consider the case A = 1. Substituting S(j, k) — C{k 1 j + 1) into (2), we get 

s-1 t-1 

Pix k = E c ( fc ^') 7V ^+E c ( fc ^' + 1 ) 7V ' ?d 

3=0 3=0 

s-1 t 

= v + ^C{k,j)N : >v + ^C{k,j)N 1 - 1 d 

3=1 3=1 

Now (5) follows by considering cases s — 1 > t, s — 1 = t or s — 1 > t. 
We now proceed to finish the proof of the theorem: 

Case 1: |A| > 1. If one of these Aj,j = 0, • • • , w — 1 is nonzero, then we see from 

(4) that — >• oo as k — >• oo; otherwise PiXk — B is a constant vector. 

Case 2: |A| = 1 and A ^ 1. If one of these Aj, j — 1, • • • ,w — 1 is nonzero, 
then we see from (4) that ||.Pi2;fc || — > oo as fc — > oo; otherwise the sequence 
PiXfc = B + A fe A , fc = 1, ■ ■ ■ , is bounded. 

Case 3: A = 1. If one of these Bj, j = 1, • • • , I is nonzero, then we see from 

(5) that ||PiXfe|| — > oo as fc — » oo; otherwise PjXfc is the constant vector u. 

Case 4: |A| < 1. We see from (4) that the sequence PiXk,k = 1, • • • , con- 
verges to the constant vector B. 

We conclude that for each i, the sequence PiXk,k — 1, • • • is either bounded 
or tends to infinity. If one of these sequences tends to infinity, then since 
\%k\ > ll-Pj^felh the sequence Xk,k = tends to infinity in norm | • | and 
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hence in norm || • ||, since the two norms are equivalent. Otherwise all sequences 
PiXk, fc = 1, ■ ■ ■ are bounded for alii = 1, • • • , m, and from which it follows that 
Xk, k = 1, • • • is bounded. This completes the proof. Q.E.D. 

The following example shows that in infinite dimensional spaces, Theorem 
1 is false. 

Example 1 For nonnegative integers n, let c„ = \n(n + 1). Define A, = i 
for c 2n < i < c 2n+1 - 1, n = 0, 1, • • • and Aj = 2 for c 2 „_i < i < c 2 „ - 1, n = 
1,2, Define a linear operator A : l 2 — > / 2 such that Ae^ = A,e i+ i,i = 
0, 1, • • •. Then the sequence of iterates A k eo = AoAi • • • Xk-iet, k = 1, 2, • • • con- 
tains subsequences that converge to 0, (e.g. A C2ri - 1 e — 1/2", n = 1,2, •••), 
subsequences that approach infinity (e.g.A C2 ™eo = 2™,n = 1,2, •••), and in- 
finitely many bounded nonconvergent subsequences (e.g. A n eo = e n for n = 
2, 4, 8, 12, 18, 24, 32, 40, • • •)• Clearly A is bounded but not compact. 

Remark 1 For T linear, Theorem 1 appeared in a 1997 unpublished article 
"On the behavior of the iterates of a matrix" of the author. 

In the case that the mapping T in Theorem 1 is linear, we can actually say 
more: 

Theorem 2 Let A : C" — > C" &e a linear map. Let \\ ■ \\ be a norm on C". 
Then for any vector x e C™, either lim^oo T k x = or lim^oo ||T fe x|| = oo or 
F < ||T fc x|| < G for sufficiently large k's, where F,G are positive numbers with 
F < G. 

Proof. 

If c = in Theorem 1, then the constant vectors B and d in its proof are 0. It 
follows that in all cases where PiXk 7 k = 1, 2, ... is bounded and not convergent 
to 0, it is bounded away from 0. If PiXk is bounded away from for some i, 
then so is Xk = T k x since \xk\ > Otherwise PiXk —¥ converges to 

for all i and hence so does xu- Q.E.D. 



Example 2 The following simple example shows that linearity is needed in 
Theorem 2: Let T : C — >• C be the map Tx = ix + c, where c e C is nonzero 
and i = Then T 4 ™(0) = and T 4 ™ +1 (0) = c for all positive integers n, 

showing that T fc (0), fc = 1, 2, • • • is neither convergent to nor bounded away 
from zero. 

Definition 1 Let X be a Banach space and 4:X4la linear operator. We 
say A has property (P) if X is a direct sum of two closed subspaces V\, V 2 such 
that (1) each Vi is invariant under A, (2) V\ is finite dimensional, and (3) there 
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exists < r < 1 and a positive integer N such that ||^4 fe a;|| < r fc ||x|| for all 
x £ V 2 and all k > N. 

Note that by Gelfand's spectral radius theorem, condition (3) above is equivalent 
to that A, as an operator on V 2 , has spectral radius less than 1. 
It is well-known that every compact operator, or more generally, Riesz operator, 
has property (P), see e.g. [1]. 

Theorem 3 Let (X, || • ||) be a Banach space and A : X — > X a bounded operator 
having property (P). Let c be a constant vector in X , and let T(x) = Ax + c 
for x £ X. Then for any vector x £ X , either {T k x, k = 0, 1, • • •} is bounded 
or Hindoo ||T fe x| = 00. Moreover, if c = 0, then either \imk^ooT k x = 
or linifc^oo ||T fe a;|| =00 or F < ||T fe a;|| < G for sufficiently large k's, where 

< F <G < 00. 

Proof. 

Let Vi,i = 1,2, N, r be as in Definition 1. Let P i} i = 1,2 be the projections of 
X onto Vi,i = 1,2 respectively. Define a norm | • | on X as 

|x| = + \\P 2 x\\. 

1 • I is equivalent to || • || and T commutes with Pi, i — 1,2. 

For any v £ X, write Vi — PiV, i = 1,2. Then for x £ X, and k > N one has 

||T fc x 2 || = \\A k x 2 +c 2 +Ac 2 + --- + A k - 1 c 2 \\ 

< r k \\x 2 \\ + \\c 2 + ■ ■ ■ + A N - 1 c 2 \\ +(r N + ■■■+ r k - x )\\c 2 \\ 

< r k \\x 2 \\ + \\c 2 + ■■■ +A N - 1 c 2 \\ + — ||c 2 ||, 

1 — r 

showing that T k x 2 , k = 1,2, • • • is bounded. By Theorem 1, \\T k x 1 \\,k = 1,2, ■ ■ ■ 
is either bounded or approaching infinity. Hence \T k x\ = ||T fe a;i|| + ||T fe x 2 ||, k = 
1, 2, • • • is either bounded or approaching infinity. Since norms | • | and || • || are 
equivalent, the same is true for ||T fe x||, k = 1,2, ■ ■ -. 

If c = 0, then c 2 = and T k x 2 — > as k — > 00, and the last part of the theorem 
follows readily from Theorem 2. 

Corollary 1 Let {X, || • ||) be a Banach space and A : X — > X a Riesz operator. 
Let c be a constant vector in X, and let T(x) = Ax + c for x £ X . Then for 
any vector x £ X, either {T k x, k = 0, 1, • • •} is bounded or lim^oo ||T fe x|| = 00. 
Moreover, if c = 0, then either \imk->ocT k x = or linifc^oo ||T fe x|| = 00 or 
F < ||T fe a;|| < G for sufficiently large k 's, where < F < G < 00. 

Corollary 2 Let (X,\\ ■ ||) be a Banach space and A : X — > X a compact 
operator. Let c be a constant vector in X , and let T(x) = Ax+c for x £ X . Then 
for any vector x £ X , either {T k x, k = 0, 1, • • •} is bounded or lim^oo ||T fe a;|| = 
00. Moreover, if c — 0, then either lim / t_ ! . 00 T k x = or lim/ £ _ ! . 00 ||T ,fe a;|| = 00 or 
F < ||T fe a;|| < G for sufficiently large k 's, where < F < G < 00. 
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Let us consider now the behavior of the sequence of averages of T: 

Ave k T(x) = hx + Tx H h T k ~ 1 x),k = 1, 2, ■ ■ ■ 

fc 

We have 

Ave fc T(x) = I (a; + Ax + ■ ■ ■ + A k ~ 1 x) + l)c + (k - 2) Ac +■■■ + A k - 2 c) 

As in the proof of Theorem 1, we may assume that A = XI + N, and that s, t 
are defined as in there. By replacing c by x, x by 0, and hence s by and t by 
s in the proof of Theorem 1, we obtain: 
For A ^ 1, 

+ Ax + ■ ■ ■ + A^x) = \ B + g A*"' 



where 



a-l ^ 

A,- = - V ——N l x, j = 0, • • • , s - 1 



(1-A)* 

B = Y r^Wx. 



s-1 



(1 - AV 

3=0 V ; 



For A = 1, 



1 1 s 

-(x + Ax + ■ ■ ■ + A k ~ 1 x) = -raitjW'- 1 !. 

By expanding A 3 we have 



(fc - l)c + (ft - 2)Ac + • • • + A fc - 2 c = T c + Ti Nc +■■■+ T^N^c 
where 

fc-2 

T 3 = J2(k - * - l)C(i, 3)^,3 = 0, ■ ■ ■ ,t - 1. 



Assume that A ^ 1. By subtracting ATo from T and using the geometric series 
formula, one gets 

k l-\ k 
~ ~ A ~ (1-A) 2 
Using the relation (i — j)C(i,j) = (j + l)C(i,j + 1), we see that 

T 1 d T 

j+1 ~ JTTrfA j 

We shall prove in the Appendix that 

T fc .7 + 1 , C(k,3)\ k -i 
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for j = 0, • • • , t — 1, where 

D(k,j, A) = ^^-[B (fc,. ? )A J + i?i(fc,j)A-''- 1 + • • • + Bj-ifafiX + 1] 



and 



Bi{k,j) = {-iy- l C{j, i) c \ k _:^y i = 0, • • • , 3 



Note that 

B i (fc,j)-^(-l) J '- i C7(j,i) 

so that D(k,j, A) approaches 1 as fc — > oo. Substituting the formulae we obtain 
thus far into AvefcT(x), with w = max{s,i}, we get for A ^ 1, 

Ave fc T(a;) = E + ^-F + G(k, A) (7) 

where 



3=0 

and e(r) = 1 for r > and e(r) = for r < 0. 
If A = 1, then Tj is equal to 

S Jtk = C(k - 2,j) + 2C(k - 3, j) + • • • + (k - j - l)C(j,j) 

We have 

Sj,k - S jtk -i = C(k - 2,j) + C(k - 3, j) + ■ ■ ■ + C(j,j) = C(k - l,j + 1) 

where the last equality follows from repeatedly applying the identity 
C(j, i + + C(j, i) = C(j + + 1). From this it follows that 

S jlk = % fc -i+C(fc-l,.? + l) 

= Sj.fc-2 + C(k - 2,j + 1) + C(k - l,j + 1) 

= Sjj+2 + C(j + 2,j + l) + ---+C(k-l,j + 1) 

= C(j + l,j + 1) + C(j +2,j + l) + ---+C(k- l,j + 1) 

= C(k,j + 2) 



(8) 
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Therefore for A = 1, 



Ave fe T(x) = \[C{k,i)x + C{k,2)Nx + --- + C{k,s)N s - 1 x] (9) 

+ \ [C(k, 2)c + C(k, 3) Ac +■■■ + C{k, t + 1) A*" 1 ^ (10) 

Case 1: |A| > 1. Suppose that t > s. Then from (7) and (8), \X k C{k,t- 1) 
dominates all other coefficients. Since D(fc, t — 1, A) — > 1 as fc — > oo and N t ~ 1 c ^ 
0, we see that ||Ave/cT(x)|| — > oo as k — > oo. The same is true if s > t since 
A s _i 7^ 0. So assume that s — t. In the trivial case s = t = 0, we have 
Ave/cT(x) = for all fc. So assume that s = t > 1. Direct checking (see 
Appendix, item 1) shows that C(fc, j)C(j, i) g^zf^ = C(k,i)C(k — i — 2,j — i) 
so that C(k,j)D(k,j,X) is a polynomial in fc. Then from (8), we have 

G(M) = E^(^i)A fc -'(^ + ^^^c) 

A fe t_1 

= —J2pi{k,j,X)A j + P2 (k,j,X)N j c 
j=o 

where for fixed A, pi,P2 are polynomials in k. Since lim^oo \p{k) = oo for 
any polynomial p(k) we see that if 

t-i 

H(k, A) = ^ p x (k, j, X)Aj + P2 (k, j, X)Wc, 

3=0 

as a polynomial in k, is identically zero, then clearly from (7) we have AvefcT(x) — > 
E, otherwise ||Ave/cT(x)| — > oo. 

Case 2: |A| = 1,A ^ 1. Suppose that t > s. Since D(k,t — 1,A) — > 1 as 
fc — > oo and N t ~ 1 c ^ 0, we see that ||AvefeT(x)|| — » oo as fc — > oo if t > 3, and it 
approaches to or is bounded if t < 2. The same is true if s > t since A s _i ^ 0. 
So assume that s = t. In the trivial case s = t — 0, we have AvefcT(x) = for 
all k. So assume that s = t > 1. As in Case 1, we consider the polynomial (in 
k) 

H(k, A) = 2 C(fc,i)A-' + M^c) 

so that 

A fe 

G(fc,A) - —H(k,X) 
k 

If degree of _ff is one or less, then we see from above that G(k, A) is bounded 
and hence AvefcT(x) is bounded by (7). If degree of H is two or more, then 
\\G(k, A)| — > oo as k — > oo and hence so is || AvefcT(x)||. 
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Case 3: A = 1. We refer to (9) and (10). If s = t = 0, then c = x = and 
Avc fc T(x) =0. If t = and s = 1, then Ave fc T(a;) = Z. If t = s - 1, s > 2 and 
N { c + N i+1 x = for alH = 0, • • • , t - 1, then Ave fc T(ir) = x. In all other cases 
lim fe ||Ave fc T(x)|| = oo. 

The previous discussions yield the proof of the following: 

Theorem 4 Let T : C" — > C" &e an aj0?ne map defined by Tx = Ax + c, where 
A is an nx n complex matrix and c a constant vector in C". Let || • || be a norm 
on C™. For any vector ieC", define 

Ave k T{x) = + Tx H h T^x) 

T/ien £/ie sequence 

Ave k T(x),k = 0,1,2,- •• 
is either bounded or linife || J 4wefcT(.x)|| = oo. 
If T is linear, i.e. if c = 0, we can say more: 

Theorem 5 Let A : C" -> C" &e /mear map. Le£ || • || &e a norm on C™. for 
any vector ieC", define 

Ave k A(x) = \ {x + Ax H h ^^x) 

Tften t/ie sequence 

Ave k A(x),k = 0,1,2,- • • 

is either (i) convergent to 0, or (ii) \\va k \\Ave k A(x)\\ = oo, or ('m^) F < 
||vlwefe^4a;|| < G for sufficiently large k's, where F,G are positive numbers with 
F < G. 

Proof. 

If c = in Theorem 4, then E = 0, and by examining its proof we see that in all 
cases where T k x is bounded and not convergent to 0, it is bounded away from 0. 
(In Case 2 of the proof, if H (k, A) is of degree one, then G(k, A) is bounded away 
from 0, and if H(k, A) is of degree 0, i.e. a constant vector, then G(k, A) — > 0.) 
Q.E.D. 

The following example shows that Theorem 5 is false if the map is not linear. 
Example 3 Let 
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and let c = | | Define T : C 2 -> C 2 by Tx = Ax + c. Let v = ^ 
Then 

a , * 1 l-i k l-i k ~ x 
Av*T(„) = - — « + ^— " 

which does not converge to 0, and has a subsequence converging to 0. 

Theorems 4 and 5 are also valid in Banach spaces. For proof, one only has to use 
these theorems and note that the average of the bounded (resp. null convergent) 
sequence T k x 2 ,k = 0, 1, ■ ■ ■ in the proof of Theorem 3 is also bounded (resp. 
null convergent). Thus we have 

Theorem 6 Let X be a Banach space. Let T : X — > X be an affine map defined 
by Tx = Ax + c, where A is an operator with property (P) and c a constant 
vector in X . Let || • | be the norm on X. For any vector x G X, define 

Ave k T(x) = \{x + Tx + ---+ T k ~ 1 x) 
k 

Then the sequence 

Ave k T(x),k = 0,1,2,- ■■ 

is either bounded or lim/j ||^?;efeT(x)|| = oo. Moreover, if c = 0, then the se- 
quence 

Ave k A(x), k = 0,1,2, • • • 

is either (i) convergent to 0, or (ii) Ymt k \\AvekA{x)\\ = oo, or (Hi) F < 
||^4wefcAx|| < G for sufficiently large k 's, where < F < G < oo. 

Recall that compact operators, or more generally, Riesz operators have property 
(P). So the above theorem is valid for these operators. 

Our last objective is to prove the following theorem. The result concerns 
AvefcT(x) in the case s = t. It shows in particular that H is identically if 
and only if x = E, the unique fixed point of T, i.e. Tx = x, so that in case (i) 
in the proof of Theorem 4 we actually have AvefeT(a;) = E for all k, not just 
Ave k T(x) — > E as k -> oo. 

Theorem 7 Let X be a fixed complex number, A ^ 0, 1. Let 

H(k) = £<7(fc,,-)A- J + ^W Nic ) 

be defined as previously. Let 1 < i < t — 1. The degree of H(k) is at most i — 1 
if and only if 

N i x= N { c-\ N i+1 c-\ I N t ~ 1 c (11) 

(1-A) + (1-A) 2 + + (1-A)*-* [ 1 
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H(k) is identically if and only if 

x = jihrf + jih^ Nc + ' ' ' + o~^~w Nt ~ lc (12) 

The proof will follow from the following discussions. 
Consider D(k,j, A) — 1. We have 

and for k > j, 

C(k - j, 2) 



C(k -j + i, 2) 

-2k + 2j - i + 1 
\k-j + i)(k-j + i- l) 

H^+ 2i ".: +1 Vi-^V7i j '- i+lN ; 



fc k 2 J \ k J \ k 

- K- 2 +^)( i +^+^ + ---)( i+i ^ ±i +^ + 

y, -2(j - j)"- 1 - (£- l)((j - z)"- 1 - (j - j + I)"' 1 ) 
m— 1 

The coefficient of in the sum above is obtained from the following calcula- 
tions: 

m — 1 m—2 

-2 J2 U - i) P U - * + ^) m - 1 ~ P + (2.7 - * + 1) E V - *) P (i - 1 + l ) m ~ 2 ~ P 

m — 2 m — 2 

= -2(j - z)™- 1 - 2( 3 - i + 1) E V - *) P C? - 1 + l ) m ~ 2 ~ P + ( 2 -7 - * + !) E ^ " ~ * + 

m-2 

= -2(j - l y m - 1 + (* - 1) e - *) p c? - * + 

p=0 

= -2(j - z)™- 1 - (i - l)((j - z)— 1 - (j - i + l)™" 1 ) 

Our first objective is to write (j — i)" 1 ^ 1 as a combination of x = l,x\ = 
i — 1, £2 = (i — 2)(i — 1), • • • , x m _i = (i — m + 1) • • • (i — 1). The proof for the 
following item 1 can be found in [3] . 

1. Let f(x) = aix 1 + • • • + a\x + a be a polynomial of degree I. Let x = 
l,xi = x — 1, X2 = (x — 2)(x — 1), • • • ,xi = (x — I) ■ ■ ■ (x — l).Then 

f(x) = diXi-\ h d x 
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where 



do = /(I) 

di = /(2)-/(l) 

d2 = ^[/(3)-2/(2) + /(l)] 

<*3 - ^[/(4)- 3/(3) + 3/(2) -/(!)] 



4 = J[/(i + l)-C(U)/(0 + " - + (-i)7(i)] 
= a; 

Moreover, for any integer n with n > I, 

f(n + 1) - C(n, l)/(n) + C(n, 2)/(n - 1) + • • • + (-l) n /(l) = 

2. Applying above to the polynomial f(i) = (j — i) m_1 , we get 
(j - i)" 1 - 1 = c m -xx m -x -\ h c x 

where 

c = u-ir- 1 

ci = (j - 2)™- 1 - (j - l)™" 1 

1 

2! 



c 2 = ^ [(j - 3)™" 1 - 2(j - 2)™" 1 + (j - I)™" 1 ] 



1 



c 3 = ^[(j-4r- i -3(j-3r- 1 +3(.7-2r- 1 -(.7-ir- 1 ] 

c m -i = -^ry [(j - m)™- 1 - C{m - 1, l)(j - m + I)"" 1 + • • • + (-l)™"^ - l)" 1 " 1 ] 
(m — lj! J 

= (-I)" 1 " 1 
Moreover, for any n > m — 1, 

(j-n-ir^-CV, l)(j-n) m - 1 +C(n, 2)(j-n+l) m - 1 +- • = 

3. Let Xi, Ci, i = 0, • • • , m — 1 be defined above. Then 

p(i) = (i-1) [(j - i) m_1 - - i + l) m_1 ] = cix 1 +2c 2 x 2 +- • -+(m-l)c m _ia; m _i 
so that 

-2(j - i)™" 1 - p(i) = -2c a; - 3c 1 x 1 (m + \)c m _ r x m -\ 

Proof. 

Note that the degree of is to — 1 and the coefficient of i m ^ 1 is 

13 



(— l) m 1 (m — 1). Applying item 1 to p(i) and using the identity C(k — 
l,i) + C(k -l,i + l) = C(k, i + 1) we get 

do = p(l) = 

di = P (2)-P(i) = (i-2) ro - l -(i-ir- l = c 1 

d2 = ^ b(3) - 2p(2) + p(l)] = (j - 3)™- 1 - (j - 2)"- 1 - [(j - 2)— 1 - (j - l)™" 1 ] = 2c 2 



^ [p(fc + 1) - C(fc, l)p(k) + ■■■ + (-iyC(k, i)p(k -« + !) + ••• + (-l) k - 1 kp(2)- 



(-l) fe p(l)] 

1 

fc! 



k [(j - k - l) m - 1 - C(k, - fc)™" 1 + • • • + (-l) k (j - l)" 1 - 1 ] = kc k 



dm-i = (-l) m - 1 (m-l) = (m-l)c m _ 1 

4. Write io = = i(i — 1), • • • , i v = i(i — 1) • • • (i — p). Note that i p is i 
times x p in item 2. This change is necessary because there is an i that was 

factored out of the summation sign in our expansion of c ^^~^f ^ — 1. 

By item 3, the coefficient of -g^ in the expansion of J^-j+f 2) — 1 is 

— Yl™=o 1 (p + 2)c p ip. Therefore the coefficient L(m,j, A) of in the ex- 
pansion of D(k,j, A) — 1 is given by 

/m—1 \ 



i=0 \ p=0 / 

m—1 j 

p=0 i=0 
™-l jp+1 

S (p+ 2)^^(1 -xy 

m — 1 

= tt-xp £ (-!) p (p + 2 ) c p ap+1 .?(j -!)■■■ (i - p)(i - A)- 7 ' -35-1 

m-l / X \ p+1 

= (-l)'(p + 2)c p j(j - 1) • • • (j - p) ( T ^-\ 

p=0 ^ ' 

= f;(-i) i+i (i+i)! Ci _iC(i,i)( T ^) 1 

2—1 ^ ' 

To emphasize that c^_i depends on j, m, we write a as P(i,j, m), so 
P(i,j,m) = i[0--i-ir- 1 -C(m-l,l)0--i) m - 1 +"-+(-ir- 1 (7-ir" 1 ] 



1 


(1 


- A)j 




1 


(1 


- A)j 




1 
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and 



/ \ 

L(m,j,A)=^(-l) i+1 (i + l)!P(i-l,j,m)C(i,i) f — 

i=l ^ 



WehaveP(0,j,m) = c = (j-l)" 1-1 , P(m-1, j, m) = c m _i = P(i, j, m) 

Cj = for i > to. We define P(i,j,m) — for negative i. We have the 
following recursive relations, which will not be used in the sequel, 

P(i,j,m) = (j-i-l)P(i,j,m-l)-P(i-l,j,m-l),m> 2. 

Also we define 

m 

M(m,j,x) = + l)lP(i - l,3,m)C(3,i)x l 

i=l 

so that 

j, A) = M(m,j, j—j) 



M(m,j,x) - ffiO'i,- •• ,i TO _i)M(m- l,.?,x) H 

+ (-im(ji, • • • , j m _i)M(m - i, + • • • + (-lr-^-iOi, • • -,j m -i)M(l,j,x) 
= joji ■ ■ ■ jm-i{m + l)x m 

where ji = j — i,i = 0, • • • , to — 1 and Hi ■ ■ ■ , j m _i) is the unsigned 
coefficient of y % in the expansion of the polynomial 

{y - ji){y - h) ■ ■ ■ (y - jm-i) 

i.e. H\ = ji H \-jm-i,H 2 = J2i<kJdk, ■ ■ ■ ,H m -i = ji ■ ■ ■ j m -i- 

In particular, we have 

M(l,j,x) = 2jx 

M(2,j,x) = j(j-l)3x 2 + (j-l)M(l,j,x)=j(j-l)(3x 2 + 2x) 
M(3,j,x) = j(j - l)(j - 2)Ax 3 + {ji+h)M(2,j,x)-j 1 j 2 M(l,j,x) 
= J(j - l)((j - 2)4x 3 + (2j - 3)3x 2 + 2(j - l)a;) 

Proof. 

For fixed m,j, 

p(x) = M(m,j,x) — H\{ji, ■ ■ ■ ,j m -i)M(m — l,j,x) H 

+ (-im(ji, • • • ,j m _i)M(m -i,j, x ) + ... + (-l^H^in, ■ ■ ■ ,j m _i)M(l,j, 

is a polynomial in x with p(0) = 0. Fix 1 < k < to — 1. We shall show 
that the coefficient of x fe in is 0. Now the coefficient is 
{-l) k+l {k + l)\C{j,k)a k , where 

a k = P(k - to) - HtP(k - l,j,m - 1) + • • • + (-l) J PT J P(fc - 1, j, m - i) + • • • 
+ (-lr^H^Pik-lJ,!) 
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Since 



y m-l_ Hiy m-2 + . . . + ( _ 1)iff . ym -i-l + . . . + ( _ 1) m-l jffm _ 1 = . . . fo-j^) 

we have for each 1 < z < k, 

jrn-X _ HljT -2 + . . . + { _iy Hi] ^-l +... + (_!) — = 

From the definition of P(i,j,m), j™ ,j™ , • • • are the corresponding 
terms (with the same coefficient) in P(fc — 1, j, m),P(k — l,j,m — 1), • • •. 
It follows that cifc = 0. 

For k = m, we have = P(m — l,j,m) = (— l)" 1 ^ 1 since P(i,j,m) = 0, 
fori > to. Thus the coefficient of x m is (-l) m+1 (-l) m - 1 (m+l)!C(j,m) = 
(to + l)joii ' ' ' jm-i- This completes the proof. 

6. Let toi, TO2, • • • , be variables. Let k,j, be positive integers. Define R{k,j) 
to be the sum of 

m i m 2 • • • m% 

where a, b, ■ ■ ■ , c are nonnegative integers such that < a, b, ■ ■ ■ , c < j and 
a + b+--- + c = j. Define R(k, 0) = 1. Define S(k, j) to be the sum of 

to"to2 ■ • • m% 

where a, b, ■ ■ ■ , c are nonnegative integers such that < a, b, ■ ■ ■ , c < 1 and 

a + b-\ Vc = j. Define S(k, 0) = 1 and S(0, 0) = 1. Then we have the 

following identity: 

^2(-lYR(p + i,q- i)S( P + i - = 

i=0 

for any positive integers p, q. 
Proof. 

Consider a typical term t = m1m\ ■ ■ ■ rrip +i resulting from the summand 
(— l) l R(p + i,q — i)S(p + i — sign disregarded. Let us write c(l) = 
a, c(2) = b, ■ ■ ■ , c(p + i) = k. We may assume that c(p + i) > 1, otherwise 
the term belongs to a previous summand; this assumption also implies 
that t does not belong to any previous (smaller i) summand. Let 

A = {j: l<j<p + i-l,c(j)^Q} 

Then the cardinality \A\ of A must be greater than or equal to i because of 
S(p+i — l, i). Also the term t appears in the expansion of the summand for 
exactly C(|A|,i) times. Next let us consider how many times t appears in 
the expansion of the next summand (— iy +1 R(p+i+l, q—i—l)S(p+i, 
t can result from multiplying a term in R(p + i + l,q — i — 1) with a term 
s in S(p + + 1). Denote the exponent of m p+ i in s by b s (p + i), which 
is either 1 or 0. 
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Denote by S* and T* the set of terms in S(p+i, i+1) and R(p+i+l, q—i—1) 
respectively. Consider the following two sets: 

Uk = {s : s e S*,b s (p + i) = k, rs = t for some r £ R*}, k = 0,1. 

Note that each s G Uk corresponds to exactly one r such that rs = t. 
Clearly U\ has exactly C(|.A|,i) elements, and these elements yields the 
same number of t's which cancel out with the previous ones because of 
the sign change. 

Each s e Uo consists of i + 1 factors from mi,- - ■ ,m p+ i-i, so Uq has 
C(|A|,i + 1) elements which yield the same number of t's. 
If C(|j4|,i + 1) = 0, i.e. \A\ < i + 1, So is empty and we are done 
since no further t's will result. If not, we consider the next summand 
(-l) l R{p + i + 2,q-i- 2)S(p + i + l,i + 2). 

Denote by S* and T-J* the set of terms in S(p + i + 1, i + 2) and R(p + i + 
2, q — i — 2) respectively. Consider the following two sets: 

14 = {s : s € Si,b s (p + i) = k,rs = t for some r G Rl}, k = 0,1. 

Note that it is b s (p + i), not b s (p + i + 1), in the above definition of 14- 
Also note that for each s e Vfe, 6 s (p + i + 1) must be or else rs = t is 
impossible. Then it is clear that |Vi| = C(\A\,i+l) and \V \ = C(\A\,i+2), 
as before. The C(|A|,i + 1) t's resulting from V\ cancel out with the 
previous the same number of t's. The same statement about Uo above 
applies to Vb and the process continues. This process will continue for at 
most p — 1 times since \A\ < i + p — 1. This proves that the t's occurring 
in J2i =0 (— l) l R{p + i,q — i)S(p + i — l,i) are all canceled out. Q.E.D. 

7. We now finish the proof of Theorem 7. Let i = t — 1. H has de- 
gree of at most i — 1 if and only if the coefficient of fc t_1 is 0. Since 
limfc^oo D(k,j, A) = 1, this amounts to 

A t _ 1 + Tr L_N^c = 

which is equivalent to 

1 — A 

So the assertion in Theorem 7 is true for i = t — 1. Suppose the assertion 
is true for some i > 1; this implies that coefficients of W in H are for 
j > i. We shall prove that it is also true for i — 1. This will complete the 
proof by induction. 
By induction hypothesis, we have 

t-i 

N 2 x = V ———Wc 

4-f (1 - a >?- i + 1 

3=1 
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By applying N to both sides t — 1 — i times, recalling that N l c = 0, we 
get 



t-i 



N 3 c 



(13) 



for all? = i, ■ ■ ■ ,t — 1. 

Substituting (13) into the formula for Ai, we find 



t-i 



j-l + 1 



-d-A) 

for all Z = i, • • • , t — 1. 

Using = jz^(Aj - iV J ' _1 a;), we get 



J-J+2 



(14) 



t-i 



E j ~ j + 1 
. (1 - AV-*+ 3 



1 



(1 - A)i 



7V J c -N l ->-x 

1 — A 



(15) 



If we write Pj for the term 



C(kJ)\-UA j + 



D(k,j, A) 



N 3 c 



1 (1-A) 2 

in H and substitute (14) and (15), we find that, for i < I < t — 1 



Pi = C{k,l)\- 1 



t-i 



- y p - l + 1 NVc+ D{k ^ x) - l N i c 

^ (I - \)p- 1 + 2 
p=i+i v ' 



(1 - Xf 



(16) 



where for / = t — 1 the sum X^-=;+i i s an empty sum, and hence its value 
is 0; and 



P i - 1 = C(k,i-l)X 



-t+1 



t-1 



4-? (i- 



^ (1 - A)J-'+ 3 1 - A 



(1 - A) 2 
(17) 



Since C(k,j) and C(k,j)D(k,j,X) are polynomials in fc of degree j, the 
terms P,-, j = 0, • • • ,i — 2 in iJ contain only fc powers of power less than 
i — 1. And since the coefficients of A; powers of power greater than i — 1 
are by induction hypothesis, we see that the coefficient of k 1 ^ 1 in H is 
equal to the limit 



lim -—(P i _ 1 +p i + ... + p t _ 1 ) 

Fix an m, i < m < t — 1. Write ^4 = A(fc, m) for the number 

C(fc,m)A- m 
k m {\ - A) 2 



(18) 
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The coefficient of the vector N m c in jr=rPj for each j, i < j < m — 1, is, 
by (16), 

which can be rewritten as 

Um-j \ 

and which by item 5 is equal to 

(-l) m -iS(m — j — l,m — j — 1)L(1, m, A)) 

where 5 is defined as in item 6, with mi = m — 1, • • ■. We also have 

-k°~ i+1 = ^- l+1 (l + R(m-j,l)k- 1 + R{m-j,2)k- 2 + --- 



(k - j) ■ ■ ■ (k - m + 1) 

+R(m - j,j - i + l)k-J +l - v ) + 0(l/k) 
j-i+i 

= R(m-j,j-i + l-z)k z + 0(l/k) 



z=0 



where R is defined as in item 6, with mi = m — 1, • • -. 
The coefficient of the vector N m c in p=rP m is 

Ak m - i+1 (D(k,m,\) - 1) 
= Ak m - l (L(l, m, A) + L(2, m, A)fc _1 + • • • + L(m - * + 1, m, A)fc-" l+l ) + 0(l/fc) 

m—i 

= A^L(m-i + l- z,m,X)k z +0(l/k) 

Note that A^ m c does not appear in Pi for Z > m; so the coefficient of JV ra c 
in the sum 

j A r (P j + ...+P t _ 1 ) 
is the same as that in the sum 

l -(P l + ... + P m ). 



k l 

Therefore the coefficient of N m c in the sum 

j A T (P i + ...+P t _ 1 ) 
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which we shall call c m is A times 

m—lm—j 

Y = ^2(-l) m ~ j ~ l+1 S(m-j-l,m-j-l)L(l,m,\) Y R{m - j, j - i + 1 - z)k* + 

j=i 1=1 z=0 



Y L(m - i + 1 - z, m, X)k z + 0(l/k) 



z=Q 

(Note that A = 0(1) so AO(l/k) = 0(l/k).) 

Now we shall show that the polynomial part of Y is a constant, i.e. the 
coefficients of k z 7 z = 1, • • • , m — i are zero. 

Fix z, 1 < z < to — i. Since k z occurs in the sum X)i=o +1 R{ m ~ 3i3 — 
i + 1 — z)k z only when j — i + 1 > z, i.e. j > z + i — 1, we see that the 
coefficient of fc z in F is L(m — i + 1 — z, m, A) plus 

m—l m— j 

Y E (-l) m ^"' +1 ^(^ - i - 1, m - j - l)L(l, m, X)R(m - j,j -i + l-z) 

j--Z + i-l 1 = 1 

m — i — z-\-l m—l 

Y L(l,m,\) Y (-l) m - j - l+1 S(m-j-l,m-j-l)R(m-j,j- i + l-z) 

1 = 1 j=z+i-l 
m — i — z-\-l m—i—z-\-l—l 

Y L(l,m,X) Y (-l) n+1 R(l + n,m-i- z + 1 - I - n)S(l + n- l,n) 

1=1 n=0 
m—i—z 

= Y L(l,m,X)-0 + L(m-i-z + l,m,X)-(-l) 
i=i 

= —L(m — i — z + l,m,X) 

where we have used item 5 for the 0. Hence the coefficient of k z is 0. 
The constant term in Y is given by the coefficient of k°, which is 
L(m — i + 1, m, A) plus 

m—l m—j 

Y E (-l) m ~^ +1 S(™ - j - 1, m - j - l)L(l, to, X)R(m - j, j-i + 1) 

j=i i=i 

m — i m—l 

= Yl m ' A ) E (-ir~ j ~ l+1 R(m - J, + l)S{m -j-l,m- j - I) 
i=i j=i 

m — i 7n — i — l 

+ n — 1, n) 

1=1 n=0 
m — i 

= Y L(l,m,X)(-l) m - l - l - 1 R(l + n,0)S(m -i,m-i-l + l) 
i=i 

m — i 

= Y (-l) m ~ i ~ l ~' 1 S(m -i,m-i-l + l)L(l, to, A) 
i=i 
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Adding the term L(m — i + l,m, A) back in we conclude that the said 
coefficient c m is 

A[L(m-i+l, m, X)-S(m-i, l)L(m-i, m, A)H \-{-l) m ~ i S{m-i, m-i)L{\, to, A)]+0(l/fc) 

which by item 5 is 

A^( m - i + 2)(— ) + 0(1/*) 

On the other hand, by (17), the coefficient of A m c in Pi_i/fc I_1 can be 
written as 

-A— (to- * + !)(— ) +0(1/*) 

Therefore the coefficient of 7V m c, i < to < i - 1, in pM-^-i H 1" p t-i) 

is 

i / \ \ m—i+l 

which approaches to 

A-' +1 
(i-l)!(l-A) m - i + 3 

as fc — )■ oo since A = A(k,m) — > A~ m /(TO,!(1 — A) 2 ) as fc — »■ oo. Since 
D(k,i — 1,A) — > 1 , c fcVlT"^ — > (i-i)! , and since iV J_1 only appears in 
Pi_i, we see that the above sentence is also valid for m = i — 1. It then 
follows from (17) and (18) that the coefficient of k 1 ^ 1 in is 

' t-i 1 1 

y - — A .. .,, jv j c -iv^ 1 ^ 

4^ (1 - AV- 4 + 3 1 - A 

Therefore 7? is of degree i — 2 or less if and only if 

t-i 

^- i .t= y - — — — N^c 

4^ (i - xy-'+ 2 

The last part of the theorem corresponds to i = 1 in the above equation. 
This completes the proof of Theorem 7. 

Appendix 



A-' +1 
(i-l)I 



We prove by induction that 



^-(l-A)^ 1 (l-A)^+ 2 + (1-A) 2 U ^^ A > 
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for j = 0, • • • , t — 1, where 

D(k,j, A) = ^-L^[B (fc, j)A'' + B.ikJ)^- 1 + ■■■ + Bj-iihJ)* + 1] 

and 

As we noted in the paper proper this is true for j = since -D(/c, 0, A) = 1. By 
item 1 below Tj can be rewritten as 

k j + 1 \ k ~iA 

(i - A)^ 1 ~ (i - xy+ 2 + (i - xy+ 2 

where 

j 

A = ^(-l)^C(fc, i)C{k -i-2,j- *)A 3 " 4 
»=o 

Taking the derivative of Tj with respect to A, we find 

d k(3 + 1) (j + + 2) A^^([(l - A)(fc - + (j + 2) A] A + (1 - A)A^) 

d\ 3 (l-A)J+2 (l-A)J+3 (l-A)J+3 

Then item 3 below proves that 



d_ k(j + l) (j + iXj + 2) C^j + ljA^ 

dA J (1-A)J+ 2 (1-A)J+ 3 U ' (1-A) 2 



This completes the induction since, as stated in the paper proper, 

1 d 



1. 



tj+1 -JTTdA Ti 



C(fc,i)C(i, = C(k, i)C(k -1-2,3-1) 



Proof. 



k\j\{k - j)\{k - i - 2)\2\ 
j\(k - - i)H\(k - j - 2)\2\(k - i)\ 

k\{k-i-2)\ 
i\(k-i)\(k-j-2)\(j-i)\ 
C(k,i)C(k-i-2,j - i) 



Q.E.D. 
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2. For k > i + 2,j > i - 1, 



(A; - i + l)C(fe, i - l)C(k - i - l,j - i + 1) + (k - i - j - 2)C{k, i)C(k -i-2,j-i) 
= (j + l)C(k,i)C(k-i-2,j-i + l). 

Proof. 



(k — i + l)C(k, i - l)C(k - i - l,j - i + 1) + (k - i - j - 2)C(k, i)C(k -i-2,j-i) 
k\(k-i-l)\ 



+(k-i-j-2) 



{k-i+ l)!(i - l)!(jfe - j - 2)!(j - i + 1)! 
jfe!(A-i-2)! 



(fc-i)!i!(fc-j-2)!(j-i)! 
ik\(k-i - 1)! 



(jfe-i)!i!(jfe- j-2)!(j-i + l)! 

- i + l)k\(k -i-2)\ 



+ ( \ k -i)\i\(k-j-2)\(3-i + l)\ 

k\(k - i - 2)\[i(k - i - 1) + (k - i - j - 2){j - i + 1)] 

(k-i)W.(k- j - 2)\(j - i + 1)\ 
fc!(fc - i - 2)\[kj + k- j 2 - 3j - 2] 

(k-i)W.(k- j - 2)\{j - i + 1)\ 
fc!(fc-»-2)!(j + l)(fc-, ? -2) 
{k - i)U\(k - j - 2)\(j - i + 1)\ 

k\(k - i - 2)\(j + 1) 
(k - i)W.(k - j - - i + iy. 
= (j + l)C(k,i)C(k-i-2,j-i + l) 

Q.E.D. 

3. Lot 

j 

A = ^(-l)^C(fc, i)C(k -i-2,j- i)A^ 4 
»=o 

and let A' = 4- A Then 

[(1 - A)(fc - j) + (j + 2)\}A + (1 - A) A A (19) 
= (.7 + l)^(-l) J+1 - 4 C(fc^)C(fc-»-2,. ? + l-*)A-''+ 1 - 1 (20) 

8=0 

Proof. We have 

A = - i)C(M)C(fc - * - 2, j - OA'"'" 1 

i=0 
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It is easy to see that 

(l-X)[(k-j)A + XA'] 

3 

= (1 - A) - i)C{k, i)C{k —i — 2,j — i)X j - 1 

i=0 

and 

3 

-X^i-iy-^k - i)C(k,i)C(k -i-2,j- + (j + 2)\A 

= -i-j- 2)C(k,i)C(k -i-2, 3 - 

i=0 

Thus (19) is the sum of 

3 

B = ^{-ly-'ik - i)C(k, i)C(k -%-2,j- 

8=0 

and 

c = -i-3- 2)C(k,i)C(k - i - 2, j - 

t=0 

The coefficient of A i+1 " t in B + C is 

{-iy- i+1 [(k-i+l)C(k,i-l)C(k-i-l,j-i+l)+(k-i-j-2)C(k,i)C(k-i-2,j-i)}, 

valid even when i = or j + 1 since by convention C(n, x) = for x < 0. 
This is equal to the coefficient of A J+1_I in (20) by item 2 above. Q.E.D. 

Acknowledgement: The author wishes to thank Professors Sing-Cheong Ong and 
Timothy Murphy for their discussions and Professors Roger Horn and Charles 
Johnson for their communications. 

References 

[1] Dowson H. R., Spectral Theory of Linear Operators, Academic Press, 1978. 

[2] Lim T.C., Nonexpansive Matrices with Applications to 
Solutions of Linear Systems by Fixed Point Iterations, 
Fixed Point Theory and Applications.Volume 2010 (2010), 
Article ID 821928, 13 pages. doi:10.115 5/2010/821928. 
http : / / www .hindawi.com/journals/fpta/2010/821928. html| 

[3] Merris R., Combinatorics, Wiley-Interscience, 2003. 

[4] Bayart F. and Matheron E., Dynamics of Linear Operators, Cambridge 
Tracts in Math., no. 179, 2009. 



24 



